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Abstract. We consider the problem of c-axis transport in double-layered cuprates, 
in particular with reference to P^S^CaC^Os+i compounds. We exploit the effect 
of the two barriers on the thermal and tunnel transport. The resulting model is able 
to describe accurately the normal state c-axis resistivity in P^S^CaC^Og+^j from 
the underdoped side up to the strongly overdoped. We extend the model, without 
introducing additional parameters, in order to allow for the decrease of the barrier 
when an external voltage bias is applied. The extended model is found to describe 
properly the c-axis resistivity for small voltage bias above the pseudogap temperature 
T* , the c-axis resistivity for large voltage bias even below T c , and the differential dl/dV 
curves taken in mesa structures. 
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1. Introduction 

The c-axis conductivity in layered superconductors is at the same time a fascinating 
puzzle and a source of interesting information on the nature of high-T c cuprates. 
Focussing on Bi 2 Sr 2 CaCu208+5 (Bi:2212), on passing from underdoping to overdoping, 
the c-axis resistivity p c changes from monotonically decreasing from T c to high 
temperature to a non-monotonic behaviour, displaying a minimum at a temperature 
usually higher than the pseudogap temperature T* [TJ IS]. This latter feature has 
led to a description in terms of tunnelling phenomena along the c-axis accompanied by 
some kind of charge localization on the ah planes [HEj- Such localization becomes even 
stronger in the so-called pseudogap phase as testified by an increase of p c and a 
change of slope in the in-plane resistivity (p a b) PQ. 

In addition to the study of the c-axis resistivity, many investigations have been devoted 
to the study of the nonlinear electrical response exhibited by the I — V characteristic 
IE], dl/dV curves [OJ EH HU H2], Pc at high applied voltages (THj or a combination 
of such measurements [H|. While in general terms it is fair to say that by means of 
such measurements one can get information both on the normal state underlying the 
pseudogap state and on the pseudogap state itself as well as on the superconducting 
state quasiparticle response, the discussion on such data is still very active [15]. 
Confining for the moment our attention to some feature of the measurements taken in 
artificially patterned mesa structures, the so called "Intrinsic Tunnelling Spectroscopy" , 
we recall some of the results that are particularly pertinent to the present work. In 
overdoped and optimally doped BSCCO mesas [HUEl], composed by several unit cells, 
the differential conductivity exhibits two distinct behaviours for voltages higher than a 
characteristic value, V > V g , and for small (vanishing) voltages, V — > 0. The (V — > 0) 
behaviour usually regarded as originating from quasiparticles, but the strangest features 
came from the large voltage behaviour indicated as "normal" resistance Rn- In fact, it 
was found that Rn was linear in T down to T — > 0, in contrast with the usual behaviour 
of p c . This feature does not seem to have been satisfactorily described. 
In two recent papers [3112] we have addressed the issue of the out-of-plane, normal state 
p c resistivities within a phenomenological model for bilayer superconductors, based on 
the existence of two energy barriers with different height and width, over which the 
electrical transport was determined by two complementary processes, namely incoherent 
tunnelling [Hj and thermal activation across the two barriers. In spite of the simplicity of 
the model, we have been able to fit a very wide set of p c data taken in BSCCO following 
the change of the temperature dependence of the c-axis resistivity with increasing doping 
5. More precisely we showed that: 

(i) For T > T* the model (which explicitly refers to the genuine normal state, with an 
electron density of states that is uniform in energy) was able to describe the p c data 
with a reduced number of fitting parameters, with a direct physical interpretation. The 
accuracy of the fit was as good in the underdoped region as in the overdoped one. 
(%%) For T c < T < T*, using the normal c-axis resistivity as determined by our model, we 
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described the loss of conductivity in the pseudogap state by a reduction of the effective 
number of charge carriers. We obtained the fraction r](T) of the "gapped" carriers that 
do not participate to the conduction along the c-axis, defined through 

* c = a c Al ~ V(T)} ^ Pc (T) = (1) 

and we found that the curves rj(T) at different dopings could be scaled onto a universal 
curve [121 HZj • 

In this paper we extend the analysis of the data by means of the proposed model, by 
extending it to nonzero voltages, thus allowing a comparison of the model with the data 
taken on the mesa structures. The main result of this paper will be a comprehensive 
and satisfactory explanation of the mesa resistivity branches, both for V — > (Rq) and 
for V >V g (R N ). 

The paper is organized as follows. In the next Section, we introduce the basic idea of 
the two-barrier model. We extend the model to finite applied voltage, and we discuss 
the limiting forms. In SecOJwe compare the model to a variety of measurements taken 
at various T on BSCCO crystals. Conclusions are reported in Sec|U 

2. Two-barrier model 




Cu0 2 Cu0 2 
< ► 

dj+d2 



Figure 1. Sketch of the energy profile used to develop the model for the small 
voltage c-axis resistivity. The different processes (thermal activation and tunnelling) 
and heights Aj and widths di of the barriers are depicted. The barriers can be thought 
of as arising from the Ca layer and the Sr-Bi-Bi-Sr block, respectively. 

As already mentioned, in two recent papers [3J HH] we have presented a 
phenomenological model for the out-of-plane resistivity p c in bilayer superconductors. 
In this model the temperature dependence of the c-axis resistivity in the genuine normal 
state (T > T*, where T* indicates the pseudogap opening) arises from the existence of 
two different energy barriers along the c-axis, taken of rectangular shape with different 
heights (Ai and A 2 , as measured from the Fermi level) and widths d\ and o^, as depicted 
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in Fig^ In this spatial energy landscape, in the normal state two mechanisms (for each 
barrier) contribute to the interlayer transport: (i) thermal activation of the carriers 
above the barrier (ii) tunnelling of the charge carriers partially inhibited by incoherent 
in-plane scattering. 

The tunnelling contribution over each of the barriers was identified as given by the model 
of Ref. jS], for which the tunnelling contribution to the c-axis resistivity is modified by the 
existence of the interplanar incoherent scattering, leading to a confinement of the charge 
carriers in the a, b planes. In this model, the in-plane phononic, inelastic scattering, 
with characteristic time r, reduces the tunnelling matrix element t c to (2r/%)t 2 c « t c 
when t is much smaller than the time between two tunnelling elements. These features 
of the model imply a relation between the in-plane and the out-of-plane resistivities. 
The overall expression for the normal state effective resistivity along the c-axis, p c>n , is 
obtained as a function of the in-plane p a b, n from the series of the two barriers [B]: 



"1 + "2 j =1)2 
i=l,2 ^1 + d 2 



t _^A _j_ g e -Ai/k B T 
apab,n + b 

-1 



(2) 

Ptun,i Pth, 

where the last equality defines the tunnelling (p tU n,i) and thermal (pth,i) resistivities 
of the i—th barrier, respectively, the term ap a b,n + b takes into account the increased 
localization due to the in-plane phononic scattering and di (i = 1,2) is the spacing 
between the layers (we have taken d\ — 3 A, &i — 12 A, the bilayer and the Sr- 
Bi-Bi-Sr block thickness, respectively). In this expression tl Qi oc exp[— 2diy / 2m*A i /h] 
(with m* = 4.6m e [TBI) an d is the height of the z-th barrier. Since the in-plane 
resistivity p a b,n is independently measured, Eq|3 contains five adjustable parameters: 
a, b, (3, Ax, A 2 , with the additional constraint that A 1 — A 2 is doping- independent. t 2 cQi 
is determined by di, A i; m* , with its prefactor absorbed in a and b. Note that a, b, (3 can 
be combined in such a way that only two of them determine the shape of the fitting 
curve, while the third affects the fit only through a scale factor. 

In order to avoid unnecessary complications, we recall that the thermal contribution to 
the higher barrier (1) was found to be negligible [3], with A 1 /k B > 2000K, so that for 
this barrier we can safely assume p\ ~ pt un ,i in the whole temperature range explored 
(T <300K). The reference expression for the effective resistivity then simplifies to: 

1 / , P2,tunP2,th , \ / q \ 

Pc,n = , , , Pl,tundl H : «2 (3) 

«1 + «2 \ P2,tun + P2,th / 

In Figure El we present a sketch of the resulting temperature dependence of the c- 
axis normal state resistivity p C) „. A thorough discussion of the parameters involved in 
equation and of their dependence on doping 5 has been given in a previous paper 
PJ. However, it is interesting to plot separately the various contributions to the overall 
resistivity, p tun<i and p thji . 

In particular, equation [3] predicts for p c n the following general behaviour: at very low 
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Figure 2. General overview of the temperature dependence of the c-axis resistance as 
calculated in the two-barrier model (black thick line), R c . n = Pen dl ^ d2 , where S is the 
transverse section. Here we have neglected (see text) the thermal conductivity across 
the higher barrier. The contributions of the first barrier (red thick line) and the thermal 
(dashed line) and tunnel (thin continuous line) contributions of the second barrier are 
separately plotted (Ri t tun — Pi,tun% f° r the tunnelling resistance of each barrier, and 
similarly for the thermal contribution). We also report the effect of increasing A 2 
(dotted lines). 



T, p2,th 3> P2,tun, so that p c>n is substantially the sum of the two tunnelling resistivities 
pijun and p2,tun, both of them linearly increasing as a function of T. As a result, p CtU 
increases linearly with T at low temperatures. At intermediate T, p 2t th < P2,tun so that 
in this temperature region the resistivity of the barrier 2 is dominated by the thermal 
activation term and is thus decreasing with T, and so is also the overall p c>n . Finally, at 
high T p 2t th — ► 0, so that p Cjn ~ p^tum again increasing linearly with T. The crossover 
from thermally activated {p c ,n ~ P2,th) to tunnel dominated behaviour {p c , n ~ Pi,tm) 
occurs for A 2 /fcsT ~1. It is interesting to note that this crossover reproduces in a 
natural way the minimum experimentally found in measurements of p c in optimal and 
overdoped samples, assigning it to a smaller A 2 with respect to underdoped samples. In 
this framework, the minimum is not a manifestation of the opening of the pseudogap, 
which occurs at a lower temperature. We stress that, as mentioned in [3], the fits of p c 
(at low voltages) only allow us to set a lower bound for Ai — A 2 . 

A second interesting feature of the model is that it naturally leads to c-axis resistance 
linear in T at low temperature, as observed in mesa structures at high voltage bias 
(that is for applied voltages V > V g , where V g is a typical voltage applied to the mesa 
structure above which the I — V curve becomes approximately linear). These data are 
reported, e.g., in Fig. 2 of Ref.[T3] and in Fig.l (inset) of Ref. |Hj- 

Within our model, as the voltage is increased the effective height of the second barrier 
decreases (see also below) while the much higher first barrier remains substantially 
unaffected. As a result, both the thermal and tunnelling resistivities of the second 
barrier become much smaller than their zero voltage values, so that the high voltage 
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(4) 



(S is the cross-section of the mesa structure) which, due to the linear dependence (as a 
function of T) of the in-plane resistivity directly explains the (so far unexplained) linear 
behaviour of Rn, as observed in mesa structure data [T3l ITIj. 

In the following, we analytically extend our model to nonvanishing voltages V, in order 
to compare the predictions with the large amount of data taken on mesa structures. 
To be quantitatively applied to measurements of the nonlinear resistance, I — V and 
dl/dV characteristics, the model requires a generalization that includes a direct effect of 
the applied voltage on the barrier height. In this case the rectangular barrier framework 
is no longer appropriate and a more accurate modelling of the barrier must be performed 
in order to take into account the decreasing of the barrier height with the external electric 
field. As depicted in FigJHJ we have chosen parabolic barriers since this choice allows 
for analytical results The energy profile along the c-axis (x direction) is then given 
by two parabolic barriers with the shape: Ui(x) = U C i — |(x — y) 2 in the interval 
< x < di, and zero otherwise. The parameter U c ^ is directly related to the previously 
defined height of the barrier A» through the relation U c> i = A« + Ep (see also figure 
EJ), while a = 8U C i/df. We stress that the choice of parabolic barriers, given the same 
widths di as in the rectangular barriers, does not introduce additional parameters. 
Now we describe the different transport processes in the finite voltage regime. We 



Figure 3. Left panel: generalized two-barrier model for the study of linear and 
nonlinear resistivity. The parabolic shape of the barrier is one of the possible choices 
that have to be selected, in order to incorporate in the model the reduction of the 
barrier height due to an external bias (finite voltage). The effect of the external 
voltage is sketched in the right panel. 



consider first the activation process. In the presence of an external electric field 8, the 

% Many other choices are clearly possible. We mention that we have performed the calculations and 
fittings also with trapezoidal barriers, obtaining the same results as for the parabolic case, with slightly 
varied numerical factors in the parameters. 
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barrier becomes U'^x) = Ui(x) — eSx. Assuming that the volume density of carrier 
n = const outside the barrier, the thermal charge current density is [T§] : 



Ji 



th,i 



eDn (l 



p i>2 e [ui(*)-E F ]/k B T dx 



(5) 



where D is the diffusion coefficient and V% is the voltage drop across the i-th barrier. 
Thus, one has: 



%i,2 



[U' r {x)-E F ]/k L 



dx 



^Ai/k B T 



%i,2 



/k B T 



dx 



(6) 



where £ is the electric field. The integral on RHS can be expressed in terms of the 
error function = - Jq e~* 2 dt, and after some straightforward calculation one gets 
the contribution of the activation to the diffusion transport current for a single barrier: 

eDn (l - e eV ^ kBT ) smh(eVi/2k B T) 



Ji 



th,i 



HGUa 
* kBT d, 



exp 



WU c ik B T 



2$ 



(AU cl -eV,) 2 -l&E F U cl 



16U ci k B T 



(7) 



In considering the tunnel resistivity, we make use of the expression 

,A l+ E F t 2 AE) 



Ji{Vi) 



[f(E)-f(E + eVi)]dE 



(8) 



e(ap ab + b) 

as for the field emission effect with the additional localization considered in reference 
jS], which we assume to be valid also in the nonlinear regime. Here t 2 ci (E) is the (energy 
dependent) transparency of the barrier and / is the Fermi function. To calculate t 2 ci (E), 
we use the same barrier's profile U(x) used to calculate the thermal contribution. For a 
generic energy level E in an external field the barrier's profile is shown in FigJSl After 
some calculation [213 EI] , one has: 



tl{E) = exp 



exp 



2V2 



%i,2 



ft 



-x z + 



ad; 





1/2 ^ 


^jx-E 


dx 



itj2m*U rJ di 



Ah 



eVj 



(9) 



for E < U ci (l — j^f-^J and t 2 ci = 1 otherwise. Here £j is the electric field on the z-th 
barrier. The total current in the series of the two barriers is the sum of the various 
contributions each with its specific values of the parameters. 

It must be noted that, in the nonlinear regime, several complications arise. First of all, 
the energy dependences of the current densities (thermal and tunnelling) is very complex 
by itself, and additionally due to the very similarity of the energy scales involved: the 
energy E, the Fermi level Ep and, in some ranges, the thermal energy k B T are all of 
the same order of magnutude. A second, more subtle complication is due to the nature 
of the double barrier: the experimentally controlled parameter is the external potential 
V, applied at the series of two barriers (assuming, as in measurements of the mesa 
structures, that one knows the number of junctions iV involved, the applied voltage 
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is just NV). However, V is not the voltage across each barrier, but the sum of the 
potential drops over the two barriers: V = V\ + V 2 . Due to the temperature and to 
nonlinear effects, the ratio V1/V2 is not constant. One can implicitly obtain each Vi by 
requiring that the total current through the two barriers be the same: J = J\ = J 2 , 
with Ji = J t h,i + Jtun,i (the cross section S is the same throughout the sample). One 
can thus solve numerically to obtain V\ and V2 as a function of V and the I(T, V) = JS 
and R C {T, V) = V/I for every temperature and arbitrary V. 

We are now in the position to perform some comparison between the model here 
developed and the data for the c-axis transport. The next Section is dedicated to 
this task. 

3. Comparison with experiments 

In this Section we report the application of the model developed in SecEl based on 
parabolic barriers, to experimental data for c-axis transport in BSCCO reported in 
the literature or measured in our laboratory. We first analyze the data for the c-axis 
resistance taken in optimally doped BSCCO single crystal mesa structures [Hj. Since 
those data were taken at small (V — > 0) and large (V = V g ) voltages, the simultaneous 
fitting of both sets of data allows a precise determination of the doping-independent 
value of Ai — A 2 . Once determined, this value will be kept fixed in the fits of all the 
other data at various dopings. 

Due to the illustrative nature of these fittings, we describe in some detail the procedure 



Figure 4. Normalized data for the c axis resistivity as reported in at high voltage, 
R c y > open circles, and small voltage, i? c ,y_»0j full dots, compared to our data for in 
a BSCCO single crystal with nominal 5 =0.255 (crosses). Continuous curves show the 
voltage bias dependence of R c , calculated on the basis of the model described in Sec|2| 
with A 2 /k B =690 K: thick black line, V ~ 0; thick red lines: V =175, 200, 220 mV; 
thin lines, from top to bottom: 45, 90, 160 mV. 



that we followed in order to reduce the number of fitting parameters. For the small 
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voltage resistivity, we have to refer to Eq|21 which contains the in-plane resistivity p a b. 
Since those data were not available in ^3], we compared the out-of-plane resistance 
with our data taken in BSCCO single crystals at various doping 5 by an eight terminal 
contact method, which allows us to extract simultaneously p c and p ab 0- As reported 
in FiglH our data for a crystal with nominal doping 5 =0.255 coincides with the R c y^o 
data in [TJ] within an overall scale factor. We then decided to use our data for p ab 
in the fitting. By adjusting A 2 , (3, a and b it is possible to fit well the normal state 
R Cj v-^o down to a temperature T* where the pseudogap opens (we come back to this 
point below, see also [3 EU E|)- For the data of R c ,v g one has to refer to the complex 
procedure decribed in SecEJ calculating the total J with the dynamic voltage partition 
between the two barriers. In this case the fit is sensitive to A x . In Eq.s I5I6I7I8I we 
have used Ep/ks =1000 K [IB]- By simultaneously fitting R Ct v_> and R c ,v g a strong 
constraint is put on the parameters. We found that it was indeed possible to fit both 
R c y^o and R c ,v g (a nontrivial result) with the choice (Ai — A 2 )/k B =2000 K, of the 
same order of magnitude but smaller than the value ~5000 K estimated in our previous 
works. We remark that no additional or hidden parameters are required, besides those 
above mentioned. In Fig. 0] we report the evolution of the nonlinear R c at various 
voltages. It is seen that the fitting of R c ,v g is obtained with V=180 mV, as compared 
to V g ~150 mV reported in [T3] . 

As a second illustration of the applicability of the model, we present in Fig. E]the fits 
of the small voltage resistivity, p c y^ , as measured by us in BSCCO single crystals. In 
these fits we have taken (Ai — A2) / k B =2000 K, as indicated by the simultaneous fittings 
above performed. Thus, once a simultaneous fit of low and high voltage resistance has 
been performed at a single doping, the fits at all the other dopings contain only four 
adjustable parameters (one of which acts as a scale factor only, see SecJ2}. As can be 
seen, the fits capture the minimum of p c . The behaviour of the parameters as a function 
of the doping has been already discussed [3]. Here, we recall that the fits depart from 
the data when the pseudogap opens. In this case from Eq. [T] one can evaluate the 
fraction rj of the carriers that no more participate to the conduction below T*. It is 
not the purpose of this paper to discuss the behaviour of 77, that has been partially 
discussed before [IB] and elsewhere |17j . However, for completeness we report some of 
the resulting rj in the inset of Fig. 0D 

As a third and final illustration we calculate the differential dl / dV curves by means of 
the nonlinear equations I5|fi|7|8l and we compare the result with the data measured 
in BSCCO mesa structures in [11] . The discussion of dl/dV data is of particular 
interest, since the "dip and hump" structure has been the subject of many discussions 
[3 HU [221 123 121] • In particular, the survival of this structure at high temperatures is 
not yet unanimously attributed to some specific mechanism, and it is often interpreted 
as a persistence of the pseudogap at high temperatures [2S] • 

Calculations of the dl JdVvs.V curves using our model are reported in Fig. [U] In 
order to compare the calculations with the data in [H], the curves have been shifted 
vertically. It is clearly seen that, within a single scale factor of order 1, our calculation 
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Figure 5. Typical fits of the c axis resistivity at small voltage with the proposed 
model, continuous lines, compared to our data in BSCCO single crystals at various 
dopings (symbols). The data for the most overdoped sample are reported in the lower 
panel with a much enlarged scale to show the minimum in p c . In the inset we report the 
fraction of charge carriers that do not participate to the conduction after the pseudogap 
opening. 

captures reasonably well the experimental behaviour, at all temperatures. We stress 
that, once a single parameter set has been chosen for one of the curves, the others are 
calculated without further adjustments. The agreement is fairly good, especially taking 
into account the very simple (albeit analytically cumbersome) model. 

4. Conclusions 

In this paper we have reviewed the two-barrier model for the electric transport along the 
c-axis in the double-layered superconductor BSCCO. We have extended the model to 
nonvanishing voltage bias, in order to compare the calculation with the measurements 
of the c-axis resistance taken in mesa structures. These measurements show rich 
and peculiar behaviour, namely the minimum in p c {T) at some doping, the crossover 
from concave upward R c y^ (T) to linear R c y g at high voltage bias, the shape of the 
dl/dVvs.V curves above T c . Nevertheless, quantitative fits to all these behaviours 
may be obtained using a simple extension of our earlier model. It seems that double- 
layered tunnelling might be responsible or co-responsible for many of the experimental 
observations. 
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Figure 6. Differential conductance dl / dVvs.V in a sample with 5 =0.25; symbols: 
experimental data after continuos lines: fits by means of the model described in 
Sec|21 The calculated curves have been scaled by a single overall factor. The curves 
have been shifted vertically for clarity. 
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